Scattering of ^Be and ^B and the astrophysical S17 factor 



o 
o 

< 
(N 



> 

(N 
O 
00 

O 

o 

CD 
I 



■ 

X 



G. Tabacaru,^ A. Azhari,^'Q J. Brinkley/ V. Burjan^ F. Carstoiu,^ Changbo Fu/ C.A. Gaghardi,^ 
V. Kroha,2 A.M. Mukhamedzhanov/ X. Tang/'Q L. Trache/ R.E. Tribble,i and S. Zhou'' 

'Cyclotron Institute, Texas A&M University, College Station, TX 11813 
'^Institute of Nuclear Physics, Czech Academy of Science, Prague-Rez, Czech Republic 

^Institute of Physics and Nuclear Engineering H. Hulubei, Bucharest, Romania 
'^China Institute of Atomic Energy, P.O. Box 215(46), Beijing 102 413, P.R. China 

(Dated: February 5, 2008) 

Measurements of scattering of '^Be at 87 MeV on a melamine (CaNeHe) target and of ®B at 95 
MeV on C were performed. For ^Be the angular range was extended over previous measurements 
and monitoring of the intensity of the radioactive beam was improved. The measurements aUowed 
us to check and improve the optical model potentials used in the incoming and outgoing channels 
for the analysis of existing data on the proton transfer reaction ^'*N(^Be,*B)^''C. The results lead 
to an updated determination of the asymptotic normalization coefficient for the virtual decay *B — > 
■^Be + p. We find a slightly larger value, Ctoti^B) = 0.466 ± 0.047 fm~\ for the melamine target. 
This implies an astrophysical factor, 517(0) = 18.0 ± 1.8 eV-b, for the solar neutrino generating 
reaction ^Be(p,7)*B. 

PACS numbers: 24.10.Ht; 25.60.Bx; 25.60.Je; 26.65. +1 



I. INTRODUCTION 

Measurements of the energetic neutrinos produced in 
*B beta decay have played a prominent role in our new 
understanding of neutrino properties (see Q, H, IE ^ Hi 
and references therein). ^B is produced in the sun by 
the ^Be(p,7)*B reaction. A good understanding of this 
reaction rate is needed in order to calculate the expected 
neutrino flux in the standard solar model Q ■ The deter- 
mination of the astrophysical S17 factor has, therefore, 
been the subject of intense experimental and theoretical 
effort over the past decade. This work has been summa- 
rized in several recent publications i . In spite of 
these efforts, there is no clear consensus on the value of 
817(0) at the desired 5% precision. Consequently, several 
new experiments are under way or planned. 

We previously reported measurements of the asymp- 
totic normalization coefficients (ANC) for ^B using 
the proton transfer reactions ^''B(''Be,*B)''Bc 9J and 
"N(^Be,SB)i3c 0. The ANCs determine the" ampli- 
tude of the tail of the overlap integral of the ground state 
wave function of ®B onto the two-body channel ''Be -I- p. 
To find the ANCs with the i4N(7Be,8B)i3c transfer re- 
action, a distorted wave Born approximation (DWBA) 
calculation is carried out and compared to the measured 
data. The DWBA calculation needs optical model pa- 
rameters (OMP) for both the incoming ''Be-I-^'^N and 
outgoing *B-|-^^C channels. Here we report a measure- 
ment of ^Be elastic scattering on a melamine (CsNgHg) 
target, where we doubled the angular range and improved 
the monitoring of the intensity of the ^Be radioactive 
beam relative to our previous measurement |lOj | . The ex- 
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tension of the angular range was done to obtain a better 
determination of the optical potential in the incoming 
channel. We also measured the elastic scattering of a ^B 
beam on a C target with the aim of checking, for the first 
time, the OMP that were used in the outgoing channel 

8B+13C. 

Below we describe the radioactive beam production, 
the experimental setups, and the procedure for the data 
reduction. We then give results of calculations for opti- 
cal potentials in both the incoming and outgoing chan- 
nels based on a double folding procedure with an effec- 
tive nucleon-nucleon interaction. We discuss the con- 
sequences of the improved secondary beam normaliza- 
tion, and compare revised DWBA calculations to the 
i''N('^Be,*^B)i3c (ja,ta from Ref. _10] in order to extract 
a new value for the ANCs and the corresponding astro- 
physical factor S17. 



II. THE EXPERIMENTS 

The ^Be radioactive beam was produced and sepa- 
rated using the Momentum Achromatic Recoil Spectrom- 
eter (MARS) [13. The primary beam was ^Li at 18.6 
MeV/ A delivered by the K500 superconducting cyclotron 
at Texas A&M University. It bombarded a liquid nitro- 
gen cooled H2 gas target at a pressure of 2 atm, with 
entrance and exit windows of 12 /xm thick Havar. A sec- 
ondary beam of ^Be at 12.5 MeV/A was filtered from 
other reaction products by MARS. The characteristics 
of the beam spot, which were measured with a 900 /im 
thick two-dimensional Position Sensitive Silicon Detec- 
tor (PSSD) placed at the MARS focal plane (the tar- 
get detector), were a spot size of 2.5 mm x 3.6 mm 
FWHM (horizontal x vertical) and an angular spread 
of 1.8° X 0.6°. The purity of the ^Be beam was 99% 
at an average rate of ~ 80 kHz. Alpha particles were 
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the primary contaminant. For a detailed description of 
radioactive beam production with MARS, see Ref. [l^ . 
Following beam tuning, the secondary target, a melamine 
foil with a thickness of 1.5 mg/cm^, was moved into the 
beam spot. Four 5x5 cm^ PSSDs were placed symmetri- 
cally around the target on an aluminum plate, as shown 
in Fig. n Detectors 1 and 2 (110 /xm thick) covered a 
laboratory angular range from 4° to 19°, and detectors 3 
and 4 (65 /xm thick) covered 16° to 30°. All four PSSDs 
were backed by 500 /im thick silicon detectors provid- 
ing particle identification spectra {AE,E). Each PSSD 
was position calibrated using a mask with 6 slots that 
were 0.8 mm wide and were spaced 8 mm apart. The 
detectors were cooled to approximately — 10°C with two 
electric thermocoolers fixed on the aluminum plate in or- 
der to decrease the inverse current in the detectors and 
minimize their noise. The assembly was placed on a XYZ 
optical table for precise positioning. 

In our previous experiments with a ^Be beam [ol llOj . 
the number of secondary beam particles was determined 
indirectly by measuring the intensity of the ^Li primary 
beam in a Faraday cup, and normalizing the yield at 
low primary beam intensities by counting the ''Be with 
the target detector. Periodically (typically once a day) 
the calibration procedure was repeated to check for any 
rate variations due to drifts in MARS power supplies. 
The primary ^Li beam intensity was substantially higher 
for the experiment on the melamine target ^3 than the 
experiment on the "'^"B target. 

Following these two measurements, we modified the ex- 
perimental setup by adding a monitor detector to count 
the radioactive beam particles directly. The beam mon- 
itoring system, which is shown in Fig. ^ used a wire 
mesh screen to reduce the secondary beam intensity and a 
plastic scintillator coupled with a photomultiplier tube to 
count the radioactive beam particles that passed through 
the target. In parallel, we ran the old monitoring system 
with a Faraday cup for the primary beam, and compared 
the results. In subsequent measurements with a high- 
intensity ^^B primary beam, we observed a difference be- 
tween the two normalization procedures. Beam heating 
reduced the density of the gas in the production target, 
causing a drop in the isotope production rate per nA of 
primary beam current and increasing the central beam 
energy. 

By scaling the heat deposition of the beam in the gas 
target, we concluded that this effect may have produced 
a small but non-negligible shift in the beam normal- 
ization during the previous ^^N(^Be,*B)^'^C experiment. 
The effect was to over estimate the number of secondary 
beam particles and hence reduce the cross section and 
the resulting ANC. In contrast, the primary beam in- 
tensity for the ^°B(^Be,^B)^Be measurement was suffi- 
ciently low to have a negligible effect. For the present 
experiment, the monitor detector was a NE102A plastic 
scintillator coupled to a photomultiplier tube. To mini- 
mize rate-dependent effects in the photomultiplier tube, 
two screens with a transparency of 9% each were added 



to attenuate the beam intensity. The yield in the monitor 
detector was calibrated using the procedure described in 
Ref. [13. 

For the ^B elastic scattering measurement, the radioac- 
tive beam was produced via the ^II(-'^°B,^B)'^H reaction 
using a 27 MeV/A ^"B primary beam on the same LN2- 
cooled gas cell. The cell contained H2 gas at 3 atm 
pressure, corresponding to a target thickness of w 10.8 
mg/cm^. Entrance and exit windows were made of 50 
/im (42 mg/cm^) Havar. A 137 mg/cm^ Al degrader was 
placed behind the gas cell to reduce the secondary beam 
energy. A 95 MeV radioactive beam of *B was focused 
at the end of MARS with a rate of about 5 kHz. The 
beam purity was better than 95%, with a particles being 
the primary contaminant. The full- width energy spread 
was limited to 1.6 MeV using momentum defining slits. 
Beam emittance was optimized using a pair of slits af- 
ter the last quadrupole in MARS. The plastic scintillator 
behind the target was used for direct counting of the 
secondary beam particles, in this case without any wire 
mesh screen. Two telescopes, each consisting of a 110 fj,m 
thick PSSD backed by a 500 /im thick Si detector, ob- 
served the secondary reaction products. The telescopes 
covered the angular range 9iab = 4° — 19°. A 1.9 mg/cm^ 
C target was used for the elastic scattering measurement. 

In both experiments, target properties such as thick- 
ness and uniformity were verified using the radioactive 
^Be and ^B beams directly by detecting beam particles 
at 0° with and without the target. The resultant energy 
loss measurements were compared to calculations with 
the computer code SRIM to extract the thicknesses. 

III. RESULTS 

A. 'Be elastic scattering 

There were two motivations for measuring elastic scat- 
tering of ^Be from melamine. The new detector geometry 
allowed us to extend the angular region for elastics which, 
in turn, helps to define optical model parameters. Fur- 
ther, it allowed us to normalize the elastic scattering yield 
directly by counting ^Be particles after the secondary tar- 
get. 

The kinematic reconstruction of the elastically scat- 
tered reaction products was performed using the energy 
and position information from the four detector tele- 
scopes. The events selected corresponded to elastic scat- 
tering of "^Be on ^''N and ^^C since the two contributions 
could not be separated. First, we identified all the events 
with {AE, E) corresponding to ^Be, then we utilized the 
correlation of ^Be energy vs. scattering angle to select 
those that were consistent with elastic scattering off ei- 
ther ^^N or ^^C. This discriminated against scattering on 
H and inelastic scattering populating excited states in ei- 
ther ^''N or ^^C. However, it was not possible to separate 
the elastically scattered events from inelastic scattering 
leading to the first excited state in ^Be at E^x = 0.429 




FIG. 1: A three-dimensional view of the detector assembly. 



MeV. We estimated this contribution using data obtained 
from inelastic scattering of ^Li on ^■^C at 63 MeV l^j to 
the analog state at Eex = 0.477 MeV. Assuming that 
the deformation lengths are equal for the analog states, 
we calculated the inelastic cross section for ^Be on -'^^N 
and ^^C at the current energy and subtracted it from the 
data. The correction was negligible at all but the largest 
angles, where it amounted to a few percent. The solid an- 
gle calculation as a function of scattering angle was done 
using a Monte Carlo simulation that took the measured 
properties of the beam spot and the geometrical spec- 
ifications of the detector assembly as input data. The 
procedure has been described in previous publications 
(e.g. see Ref. [12]). Figure |21 shows the resulting angular 
distribution corresponding to elastically scattered ^Be on 
and 12 C. 

The angular distribution predicted from the optical 
model parameters used in Ref. T?! is compared to the 
data in Fig.|21 The Monte Carlo calculation was used to 
provide the proper angular distribution that takes into 
account the finite angular binning in the data. 

The results of our measurements are compatible with 
those reported in Ref. |0 at small angles and with the 
predictions of the optical model calculations done at that 
time (dotted curve). However, the new experimental 
data fall above the predictions at larger angles, suggest- 
ing a smaller absorption than was assumed in Ref. . In 
order to obtain a better description of the elastic scatter- 
ing, calculations were carried out with a range of new pa- 
rameters. These were also used for the entrance channel 
to generate DWBA predictions for the i''N(^Be,®B)i3c 
proton transfer reaction. 

The optical parameters used in Ref. were based 



on results from an analysis of elastic scattering of loosely 
bound p-shell nuclei 16], which demonstrated that the 
data can be described with double- folded potentials. The 
potentials quoted in Ref. |l6l ] were obtained from calcu- 
lated nuclear matter densities folded with an effective 
nucleon-nucleon interaction (JLM, |i3]), smeared (two 
range parameters, ty and tw) and renormalized (two 
strength parameters, Ny and Nw) to produce: 

UoFir) = NvVir,tv)+iNwWir,tw). (1) 

The calculations for previous ^Be studies were done using 
the JLMl effective interaction with standard range pa- 
rameters: ^^ = 1.20 fm, ivi^ = 1.75 fm, and average renor- 
malizations A''y=0.37, Nw~l-00 (for details see Ref. 
and references therein). These parameters served as the 
starting point for the new calculations. Elastic scattering 
of ^Be at 87.7 MeV on i**N and ^^C were calculated in the 
center of mass frame, then transformed into the lab frame 
and added with weights 1.0 and 0.5, respectively, equal 
to the ratio of ^''N to ^^C nuclei in the melamine. The 
resulting curve was "smoothed" using the Monte Carlo 
code described above. The parameters for the folding 
potential were varied simultaneously and identically for 
both target nuclei. This approach is supported by the 
fact that both target nuclei are well bound and have sim- 
ilar densities in the surface region and by experiments we 
have carried out with melamine and C targets using other 
radioactive beams, such as ^^N ®B (present experi- 
ment) and "F 0. The extended angular coverage of the 
present data was still not sufficient to attempt an optical 
model fit with free parameters. Rather, the two nor- 
malization and two range parameters were varied. The 
parameters for various calculations and the reduced 
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FIG. 2: (color online) Angular distribution for elastic scat- 
tering of '^Be on and ^^C. The points are experimental 
data after subtraction of the inelastic scattering contribution. 
The solid (red) line is the new calculation, and the dashed 
hne (blue) is the previous calculation from Ref. Both 
calculations are smeared to account for the finite angular res- 
olution. 



values obtained by comparison to the data are presented 
in Table n 

Four entries in the Table (A, B, C and H) were ob- 
tained by adjusting the renormalization of the real and 
imaginary parts of the potential. The smearing ranges of 
the interaction, ty and tw, were adjusted for three cases 
(D, E, F), and the density dependence was adjusted in 
one case (G) where the JLM2 interaction was used. The 
best results were obtained for cases D, E, G and H. The 
small differences between the values show that it is dif- 
ficult to choose a "best solution" . Rather, we did DWBA 
calculations for the ^^N(^Be,^B)^'^C transfer reaction for 
the four most promising potentials. The results are com- 
pared to the previous calculations in Table Q] 

A far/near decomposition of the scattering amplitudes 
shows that the observed angular range covers the region 
of Fraunhofer oscillations generated by the interference 
of the two components (see Fig. ^ . Their crossover is 
around 20°, and at larger angles the far component be- 
comes dominant. But in the region included in our mea- 
surements, the interference is still important. The cal- 
culations show that after about 60° the angular distri- 
butions develop a rainbow type pattern, typical for the 
cases found recently in our ^'^Li elastic scattering data 
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FIG. 3: (color online) Far side/near side components shown 
compared as ratios to Rutherford scattering for the potential 
H in Table U The near side components are the dashed lines, 
the far side components are the dot-dashed lines, and the 
totals are given as the solid lines. 



B. B experiment 

A similar analysis was used for the ^B elastic scattering 
on a natural C target. The resulting angular distribution 
is shown in Fig. 0] where it is compared to calculations 
made with the folded potentials using the average pa- 
rameters tv = 1.20 fm, tw — 1.75 fm, Ny = 0.37 and 
Nw = 1.00. The solid (dashed) line shows the results af- 
ter (before) smoothing with the Monte Carlo calculation. 
The ^B density used in the folding procedure was that 
calculated in jl^] using the correct ANC for the last pro- 
ton. Due to the limited angular range of the data, we did 
not attempt to produce a better fit by adjusting parame- 
ters. Based on the similar densities for ^^C and ^'^C and 
on results found in cases where scattering on both ^^C 
and "'^^C were measured, we assume that the parameters 
(ty, tw, Ny, Nyy) extracted for the natural C target are 
valid for the ^B-I-^'^C channel in DWBA calculations of 
the transfer reaction. 
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IV. IMPLICATIONS FOR i*N(^BE,®B)"C 



Based on the present analysis, the value of the 
^^N(^Be,*B)^'^C proton transfer reaction cross section 
from our previous measurement ('l(/|) was increased, 
(j%.p{new) = 1.055 cr*^p(oW), to account for the differ- 
ence of 5.5% found in the absolute normalization. Fol- 
lowing the original publication the ratio between the 
ANCs for the lpi/2 and 1^3/2 components in the wave 
function of the ground state of ^B was found to be 
J2 ^ C2(8B,pi/2)/ C2(8B,P3/2) = 0.125(20) from the 
mirror neutron transfer reaction -'^•^C(^Li,^Li)^^C [isj . 
The previous result had assumed 5^ = 0.157 based on 
theoretical calculations. The new cross section, the value 
of 5^ from experiment, and the new optical model param- 
eters have been used to find revised values of the ANCs 
for **B ^''Be + p from the ^^^{^Be,^BY'^C proton trans- 
fer reaction measurement [To|. 



FIG. 4: (color online) The angular distribution for elastic 
scattering of *B on C. The dashed line is the calculated cross 
section and the solid line is the result after accouting for the 
finite angular resolution. 



In order to extract the ANCs, four terms must be cal- 
culated to account for the possible proton transitions: 
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where the cr's are the calculated DWBA differential cross 
sections for proton transfer from the P3/2 and Pi/2 or- 
bitals in ^^N to the P3/2 and pi/2 orbitals in ®B, the 
bij's are the ANCs for the single particle orbitals used 
in the DWBA calculation, and the Cnj^ ij's and Csb Zj's 
are the ANCs for ^^^N ^^C + p and '^B ^ '^Be + p, 
respectively. The ANCs, C?4^ ^^^^ = 18.6(12) fm"! and 

C?4^ p^^_^ =0.93(14) fm-\ were measured in HHIH. The 

calculations were done with the code PTOLEMY '2i|. 
The results of the calculations are given in Table m where 
the value calculated and shown in column 7 is the quan- 
tity in the last square bracket in Eq. (01, ^dw, inte- 
grated over the angular region ^c.m. = 4°- 25° to match 
the data. This quantity contains the DWBA cross sec- 
tions weighted with the ANCs for ^^N, the single particle 
ANCs calculated for the appropriate Woods-Saxon pro- 
ton binding potentials in *B and ^^N and the mixing ratio 
in the ground state of '^B, 5^ . Since the reaction is periph- 



eral, the results do not depend on the geometry assumed 
for the proton binding potentials, which are chosen to be 
Woods-Saxon shape with depths adjusted to reproduce 
the experimental proton binding energies in ^B and ^^N, 
respectively. The results shown were calculated using the 
reduced radius ro = 1.20 fm and the diffuseness a = 0.60 
fm, and the same spin-orbit term as in Ref. [Toj |. The 
exit channel parameters were fixed to the previous val- 
ues. Calculations were done at Eiab = 83.5 MeV, the 
energy of the previous experiment with the four optical 
model sets in Table U for the entrance channel that have 
the lowest x^, B, D, E, G, and H. In column 8 we give 
the ratio of the present calculations to the same quantity 
calculated in Ref. The average of the four results, 

weighted by the x^'s, gives the ratio (7^) = 0.968 ±0.047. 
The new ANC, calculated with the relation 

Ch^p,,2^new) = (1.055/ {n))Cij,^^^^^{old), (3) 
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is Cir, (new) = 0.414 ± 0.041 fm-i. The overall un- 

certainty contains contributions from statistics (2.6%), 
absolute normalization of the cross section (5%), input 
parameters in the Monte Carlo simulation of the exper- 
iment (1.4%), and uncertainties in the ANC for the ^^N 
vertex (6.4%). The contribution of each of these factors 
remains the same as in the original analysis [lO] . The un- 
certainty due to the optical model parameters was taken 
from the standard deviation of the calculated cross sec- 
tions (column 7 or 8 in TablcQJ and is 5% compared with 
the previous value of 8.1%. 

The relation between the ANCs and the astrophysical 
factor S'i7(0), in eV-b, is H 

SMO) = 38.6 [Cl^^ + Cl^) = i^-QCl.S^ + 6^) (4) 

Using the new value of the ANC we find 5i7(0) 18.0 ± 
1.8 eV-b. This value is very close to the value obtained 
from the reaction ^°B(^Be,^B)^Be in Ref. 9] where we 
found 5*17(0) = 18.4 ±2.5 eV-b. The weighted average of 
the two results is 5*17(0) = 18.2 ± 1.7 eV-b. 

V. CONCLUSION 

Elastic scattering of ^Be at about 12 MeV/A has been 
measured over an extended angular range on a melamine 
target. The results provide a better determination of the 
optical model parameters used for the entrance channel 
of the i**N('^Be,*B)i3c reaction. For the first time elas- 
tic scattering of ^B was measured on a C target, thus 
allowing for a check the optical model parameters used 
for the exit channel in the DWBA calculation. In the 
measurement of the ^Be elastic scattering, we directly 
counted the secondary beam particles. This resulted in a 
5.5% increase of the transfer reaction cross section from 
Ref. which used an indirect method to obtain the 
secondary beam intensity. We also used the mixing ratio 
between the lpi/2 and IP3/2 components from a (^Li,^Li) 



measurement [l5j, rather than a theoretical prediction. 
These improvements lead to the revised value of the ANC 
for ®B ^ ■^Be -I- p from the "N(^Be,SB)i3c reaction 
of (7^(^5,^3/2; new) = 0.414 ± 0.041 iuT^ , resulting in 
Clti'^B-new) = C^^^^ + C^^^^ = 0.466 ± 0.047 fm"!. 
This, in turn leads to a larger value for the astrophysical 
S factor for the '^Be(p,7)8B reaction, 5*17(0) = 18.0 ± 1.8 
eV-b. This new value is very close to the one from the 
same reaction on the ^°B target 0. Averaging the two 
results, we obtain 5i7(0) — 18.2 ± 1.7 eV-b from the pro- 
ton transfer reactions. 

Our result for 5i7(0) is a bit over 2cr lower than the 
extrapolation of the most recent and precise direct mea- 
surement of ^Bc(p,7)*B by Junghans a/. 7J. Our cen- 
tral value is about 1.5cr lower than the average central 
value obtained by Cyburt et al. ^ in a recent analy- 
sis that uses all of the best available capture data, under 
the assumption that they are independent. Including the 
uncertainty quoted by Cyburt et al. our results are con- 
sistent at the 1(7 level. We do not understand the reason 
for the discrepancy between our ANC result and the ex- 
trapolated value from Junghans et al.. However, we note 
that direct measurements with both radioactive beams 
and targets and indirect measurements continue to be 
carried out on this important proton capture reaction. 
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TABLE I: The optical model parameters, the corresponding per degree of freedom for the ^Be elastic scattering fits, the 
calculated DWBA cross section for the ^''N(^Be,*B)^^C reaction (see text), and the ratio of the calculated DWBA calculation 
to that in Ref. 10] . 
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